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FOREWORD 

This  report  was  prepared  by  the  Fuel  Branch  of  the  Fuel ,  Lubrication, 
and  Hazards  Division,  Air  Force  Aero  Propulsion  Laboratory,  under  Project 
3048,  Task  304805. 

The  experimental  data  used  as  a  basis  for  this  report  is  from  tests  con- 
t  oted  by  the  Society  of  Automotive  Engineers  Technical  Committee  E-r  in 
June  1969.  Raw  data  was  reduced  to  the  final  form  presented  herein  by  a 
group  within  the  SAE  Committee  and  by  a  team  at  General  Electric  Company, 

E  vend  ale,  Ohio. 

Some  of  the  items  compared  in  this  report  were  commerical  items  that 
were  not  developed  or  manufactured  to  meet  Government  specifications  and 
were  not  necessarily  intended  for  the  service  considered  in  this  report.  Any 
failure  to  meet  the  objectives  of  this  study  is  no  reflection  on  the  value  of  these 
items  for  other  service,  nor  should  the  conclusions  of  this  report  be  construed 
as  statements  of  the  manufacturers'  abilities. 

The  analysis  described  in  this  report  was  conducted  from  September  1969 
to  February  1970  at  the  Air  Force  Aero  Propulsion  Laboratory,  Wright- 
Patterson  AFB,  Ohio  45433. 

The  author  appreciates  and  acknowledges  the  assistance  rendered  by  the 
following  outside  of  the  AF  Aero  Propulsion  Laborrtory  and  the  SAE  Com¬ 
mittee:  Mr.  C.  Fetter  of  the  Digital  Computation  Directorate,  Aeronautical 
Systems  Division,  for  guidance  in  applying  the  data  plotting  routine  "GP;" 

Mu.  Mary  Lum  of  the  Operations  Analysis  Office,  AB  Logistics  Command, 
for  reviewing  and  commenting  on  the  approach  and  on  the  analysis  criteria; 
and  Mr.  Charles  Stanforth  and  others  at  the  General  Electric  Company  for 
their  help  in  reducing  the  data  to  a  final  form  for  analysis. 

This  report  was  submitted  by  the  author  23  March  1970. 

This  technical  report  has  been  reviewed  and  is  approved. 

flaLu  t 

ARTHUR  V.  CHURCHILL 

Chief,  Fuel  Branch 

Fuel,  Lubrication,  and  Hazards  Division 

Air  Force  Aero  Propulsion  Laboratory 
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ABSTRACT 

This  report  describes  a  computerized  statistical  analysis  of  test  data  from 
engine  smoke  measurements  conducted  by  the  Society  of  Aut  >motive  Engineers 
Technical  Committee  This  Committee  was  organised  to  develop  a  reason¬ 

ably  simple,  precise,  and  universally  acceptable  standard  for  measuring  ex¬ 
haust  smoke  from  aircraft  engines.  The  analysis  indicated  that  the  Committee's 
test  data  can  be  used  to  arrive  at  statistically  meaningful  conclusions  about 
four  measuring  system  parameters.  "Whatman  No.  4"  was  found  to  be  superior 
to  "Millipore  SM"  as  a  filtering  medium  in  this  application.  All  three  reflec- 
tometers  tested  were  found  to  produce  equivalent  results.  White  reflectoneter 
background  shade  was  found  to  have  slight  superiority  over  black,  yet  black 
(i.  e. ,  absolute  reflectance  less  than  5%)  was  recommended  as  a  safeguard 
against  unknown  factors.  The  lower  sampling  flow  rate  (0,  0041  sets)  w>-s  found 
to  have  produced  slightly,  yet  consistently,  higher  smoke  density  readings  than 
the  higher  flow  rate  (0. 0085  sefs)  tested. 
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SECTION  I 
INTRODUCTION 

There  are  numerous  systems  in  use  for  measuring  exhaust  smoke  from 
aircraft  engines.  Most  consist  of  drawing  an  exhaust  smoke  sample  through  a 
filter,  measuring  the  light  reflection  from  the  resultant  spot,  and  then  com¬ 
paring  this  to  the  light  reflection  from  some  standard. 

Unfortunately,  the  many  details  of  this  seemingly  simple  procedure  have 
never  been  standardized.  Results  from  different  systems  are  not  readily  com¬ 
parable,  and  the  inherent  precision  of  most  of  these  systems  has  never  been 
defined. 

Technical  Committee  E-31  of  the  Society  of  Automotive  Engineers  (SAE) 
was  established  to  cope  with  this  problem.  Its  purpose  was  to  prepare  a 
reasonably  simple,  precise,  and  universally  acceptable  method  for  measuring 
exhaust  smoke  from  aircraft  engines. 

In  June  1969,  Committee  E-31  compv  ;nded  a  preliminary  standard  and 
conducted  tests  to  examine  the  parameters  of  this  proposed  scheme.  A  brief 
description  of  the  test  program  and  the  procedure  evaluated  are  given  in 
Appendix  I. 

"SN"  is  the  dimensionless  term  proposed  for  use  in  quantifying  smoke 
emission.  Some  of  the  test  program  raw  data  was  reduced  to  SN  by  an  analysis 
group  within  the  Committee.  Additional  data  was  later  reduced  by  a  team  at 
General  Electric  Company  (GE),  Evendalc,  Ohio.  The  data  reduction  procedures 
of  the  two  groups  differed  somewhat.  Both  are  described  in  Appendix  H. 

An  analysis  of  this  reduced  test  data  is  presented  in  this  report. 

This  analysis  of  reduced  data  was  undertaken  to  answer  the  following: 

How  much  did  each  parameter  influence  the  measurement  of  smoke  In 
comparison  to  all  other  parameters  Investigated? 

if  a  parameter  did  have  effect,  which  value  of  the  parameter  produced 
the  best  results  ? 
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SECTION  II 
ANALYSIS  PROCEDURE 

Table  I  contains  all  reduced  (SN)  data,  the  raw  material  of  this  analysis. 
Each  data  column  contains  data  for  a  combination  of  four  explicit  parameters: 

#  Filter  Medium:  Whatman  No.  4  or  Millipore  SM,  plain  white 

#  Reflectometer:  MacBeth  Model  NB-100R;  W.  W.  Welch  "Dens ichron, " 

Model  Oi.e;  or  Photovolt  Model  610 

£  Sampling  Flow  Rate:  0. 0041  standard  cubic  feet  per  Becond  (scfs) 

or  0.  0085  scfs 

%  Reflectometer  Background  Shade:  Black  or  white 
The  group  which  reduced  each  data  column  (SAE  or  GE)  is  also  noted  in  Table  I. 

The  rows  of  Table  I  are  numbered  15  through  56  in  keeping  with  the 
numbering  system  established  during  the  tests.  Each  of  these  42  rows  repre¬ 
sents  different  engine  conditions  coupled  with  values  of  parameters  other  than 
the  four  noted  above.  (Appendix  I  contains  a  complete  list  of  parameters.) 

This  is  an  important  point  that  largely  dictated  the  analysis  method  more  than 
four  parameters  were  varied  during  the  tests.  Consequently,  it  is  not  possible 
to  make  column  comparisons  unless  all  columns  being  compared  contain 
exactly  the  same  towr.  not  just  the  same  number  of  rows. 

Initially  it  seemed  possible  to  draw  statistically  valid  conclusions  about 
six  parameters.  More  detailed  scrutiny  revealed  that  this  unfortunately  was 
not  posaible.  There  was  not  enough  data  to  statistically  examine  any  param¬ 
eter  other  than  the  four  explicitly  noted  aa  column  headings  in  Table  I. 

Column  "sets'  were  established  to  overcome  the  lack  of  identical  test  con¬ 
ditions  from  row  to  row.  A  set  is  any  number  of  columns  all  of  which  contain 
the  same  rows.  Since  all  the  columns  of  Table  I  do  not  all  include  the  same 
rows,  forming  a  set  was  necessarily  a  compromise  between  getting  as  many 
points  per  column  as  possible,  while  including  as  rasay  columns  ss  possible  in 
the  set.  For  example,  see  Tables  II  through  VI. 
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TABLE  I  REDUCED  OATA  (SN)  FOR  ANALYSIS 
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The  following  four  space  code  was  devised  to  identify  the  parameter  values 
of  any  column: 

•  First  space  identifies  filter  medium 

W  -  Whatman 
M  -  Millipore 

9  Second  space  identifies  reflectometer 

M  -  MacBeth 
D  -  Densichron 
P  -  Photovolt 

•  Third  space  identifies  sampling  flow  2  +e 

4  -  0. 0041  scfs 
8  -  0. 0085  scfs 

•  Fourth  space  identifies  reflectometer  background 

B  -  black 
W  -  white 

For  example,  WM8B  is  a  column  containing  reduced  data  taken  on  Whatman 
filter  medium,  at  the  h';h  (0. 0085  scfs)  flow  rate,  with  the  resultant  spot  read 
with  the  MacBeth  meter  using  the  black  background.  Note  that  this  code  does 
not  reveal  which  or  how  many  rows  are  includeu. 

Each  of  the  four  parameters  was  analyzed  independently  of  the  other  three. 
All  the  column  sets  for  one  parameter  constitute  a  "series. "  No  one  of  the  *>ur 
ser<°s  contained  all  the  data  of  Table  I,  but  each  utilized  at  least  90%  of  that 
reduced  data 

The  computer  routine  CORRE1,  included  as  Appendix  III,  was  written 
around  two  existing  suLroutlnes  for  this  analysis.  Figure  1  is  a  typical  print¬ 
out  of  this  program  for  a  single  column  set.  The  computer  outputted  all  input 
data,  made  correlation  piots  (scatter  diagrams),  and  computed  the  following: 

9  Mean  (M)  of  each  column 

9  Standard  deviation  (30)  of  each  column 
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•  Standardized  mean  and  standardized  standard  deviation  for  each 
column 

•  Coefficient  of  variation  (CV)  of  each  column 

•  Correlation  coefficient  (r)  of  each  column  pair  specified 

These  quantities  are  defined  and  explained  in  Appendix  IV. 

There  was  no  preanalysis  attempt  to  correct  cr  exclude  suspect  data. 
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Figure  1.  Typical  Print -Out  of  Analysis  Routine  CORRE1  for  One  Set 
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SECTION  III 
RESULTS 


1.  GENERAL  CONSIDERATIONS 

The  results  of  this  analysis  are  tabulated  in  Tables  II  through  VI. 

The  number  of  data  points  per  column  is  the  most  important  statistical 
indicator  of  confidence.  About  15  points  per  column  was  generally  the  minimum 
number  that  produced  good  results. 

The  data  had  to  be  examined  closely  after  results  were  computed.  In 
several  cases,  what  appeared  to  be  poor  results  was  actually  attributable  to 
just  a  few  "odd  points,"  that  is,  deviations  from  whatever  trend  was  established 
by  the  rest  of  the  data  in  a  column  pair.  In  such  cases,  these  few  points  were 
corrected  to  a  value  that  seemed  probable,  and  the  results  were  recomputed. 
Tables  II  through  VI  contain  only  results  computed  with  uncorrected  data.  Any- 
corrected  results  are  noted  and  listed  in  the  "Comments”  column  of  each  Table. 

Comparisons  can  be  made  within  column  sets  only.  This,  as  stated  pre¬ 
viously,  is  because  of  the  lack  of  identical  test  conditions  from  set  to  set. 

Even  though  two  columns  may  bear  the  same  column  identification  codes,  they 
are  generally  not  identical  if  they  appear  in  different  sets;  the  rows  and  number 
of  rows  comprising  each  set  are  different. 

2.  ANALYSIS  CRITERIA 

Two  types  of  criteria  were  used  to  meet  the  objectives  of  this  analysis. 
Influence  criteria  were  used  to  determine  how  mucti,  if  any,  influence  each 
parameter  had  on  the  measurement  of  smoke.  Superiority  criteria  were  sub¬ 
sequently  used  to  determine  which,  if  any,  value  of  a  given  parameter  pro¬ 
duced  better  results. 

The  criteria  for  a  parameter  to  have  had  significant  influence  were: 

•  AM  (difference  in  column  means)  >  10%  of  the  lower  M 

•  ACV  (difference  in  column  coefficients  of  variation)  >  10%  of  the 

lower  CV 
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•  r  <  0. 990,  for  a  nontranslated  function 

Anyone  of  these  had  to  be  satisfied  for  a  parameter  to  be  considered  as  having 
significant  influence.  These  distinction  criteria  are  subjective.  They  were 
based  on  the  author's  preliminary  survey  of  the  computer  computations  and  on 
the  belief  of  several  Committee  E-31  members  experienced  in  smoke  measure¬ 
ment  that  the  level  of  significance  of  the  system's  results  was  about  3  SN  in  30. 

It  was  important  qualify  the  r  <  0.  990  criterion  as  being  valid  only  for 
a  "nontranslated  function.  "  The  definition  of  r  considers  dispersion  of  data,  as 
well  as  deviation  in  slope  of  the  data  regression  line  from  the  slope  of  the 
perfect  correlation  line  y  =  x.  However,  r  does  not  consider  the  effect  of  a 
translated  function  y  =  x  +  k.  Figure  2  shows  a  perfectly  translated  function. 

The  correlation  coefficient  for  both  it  and  the  perfect  correlation  line  is  1. 0. 
This  translation  phenomenon  appeared  fairly  frequently  in  the  correlation  plots. 

Two  superiority  criteria  were  used  to  distinguish  between  values  of  a  given 
parameter.  The  best  parameter  value  was  the  one  that  displayed: 

•  The  largest  M 

•  The  smallest  CV 

These  criteria  are  desirable  from  purely  mathematical  considerations  of  pre¬ 
cision.  They  are  also  desirable  criteria  considering  the  nature  of  smoke  mea¬ 
surement  and  the  definition  of  SN.  The  SN  scale  is  mathematically  defined  from 
0  to  100.  When  smoke  spots  are  rated  in  units  of  optical  density,  SN  values 
are  most  precise  at  the  scale  midpoint,  SN  -  50.  This  is  because  the  expression 
for  SN  in  terms  of  optical  density  is  a  logarithmic  function.  Also,  the  Deed  for 
precision  in  smoke  measurement  is  greatest  at  thrt  value  of  SN  corresponding 
to  the  thresh  Id  of  smoke  visibility.  Though  this  value  is  far  from  being  well 
defined,  all  work  to  date  Indicates  that  it  is  within  SN  of  *0  to  35  (References 
1  through  3).  Consequently,  the  best  vaiue  of  a  parameter  is  not  only  the  one 
that  produces  the  least  deviation  with  respect  to  the  mean  (minimizes  CV),  but 
the  one  that  tends  to  increase  the  mean  toward  SN  =  50. 
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Figure  2.  Typical  Perfectly  Translated  Correlation  Function 


3.  USE  OF  DATA  REDUCED  BY  GE  VS  SAE  GROUPS 

First  it  was  necessary  to  determine  if  the  use  of  two  data  reduction  groups 
with  somewhat  different  methods  had  .  .gnificantly  influenced  the  results.  Since 
both  groups  prepared  what  should  ideally  have  been  identical  data  w  ith  the  Photo- 
volt  reflectometer,  these  data  were  used  to  investigate  the  possibility  of  in¬ 
fluence.  Table  11  contains  the  results  of  this  comparison. 

Table  II  contains  four  column  seta,  each  with  one  column  pair.  A M  was 
insignificantly  small  with  all  four  sets,  but  ACV  appeared  to  be  significantly 
large  In  Set  2.  and  r  appeared  to  be  significantly  low  in  Sets  1  and  2. 

However,  the  correlation  ploU  reverJed  that  three  of  the  11  points  per 
column  in  Set  2  were  odd.  Set  1  also  displayed  3  odd  points  in  ite  correlation 
plot  of  22  points.  Correction  of  these  .xkl  points  made  both  suspect  r  greater 
than  0.  990.  and  reduced  the  ACV  of  Set  2  to  below  4.  2. 
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It  was  concluded  that  no  significant  difference  existed  between  the  Photo¬ 
volt  data  reduced  b  >  GE  and  that  reduced  by  SAE.  For  the  purpose  of  this 
analysis,  this  conclusion  was  taken  as  general  proof  of  the  identity  of  results 
from  the  two  data  reduction  methods. 

The  GE  and  the  SAE  Photovolt  data  were  not  mixed.  The  GE  Photovolt  data 
was  chosen  for  the  remainder  of  the  analysis  simply  because  that  team  had 
produced  more  t.  jints.  The  analysis  was  then  based  on  all  data  in  Table  I  except 
the  SAE  Photovolt  data. 

4.  EFFECT  OF  SAMPLING  FLOW  RATE 

Results  from  the  five  column  sets  that  constituted  the  sampling  flow-rate 
influence  series  are  tabulated  in  Table  III. 

The  Am  and  Acv  of  all  10  column  pairs  were  insignificantly  small,  but 
the  r  value  of  seven  of  these  pairs  was  less  than  the  0. 990  criterion.  There 
was  one  odd  point  in  one  of  these  seven  pairs,  but  even  after  correction  the  r 
value  was  still  significantly  low.  The  plots  of  two  of  the  other  three  pairs 
demonstrated  slight  correlation  function  translation,  indicating  that  their  r 
values  are  deceptively  high. 

This  small  but  significant  and  reasonably  consistent  lack  of  correlation 
indicated  that  the  sampling  flow  rate  had  a  small,  but  significant,  influence 
on  smoke  measurement. 

There  Is  also  a  consistent  trend  in  the  AM  column  of  Table  III.  The 
higher  flow  rate  (0. 0085  eefs)  produced  lower  average  SN  with  all  10  pairs 
by  0. 9  to  3.  3. 

None  of  the  10  Acv's  are  significantly  large,  so  neither  flow  rate 
appeared  t<~  have  intrinsic  superiority. 

5.  EFFECT  OF  REFLECTOMETER  CHOICE 

The  results  of  this  series  for  the  parameter  values,  MacBeth,  Densichron, 
and  Photovolt  reflectomet  "8,  are  given  in  Table  IV.  The  four  sets  of  this 
series  consisted  of  18  column  pairs. 
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TABLE  IZ .  RESULTS  OF  REFLECTQMETER  INFLUENCE  INVESTIGATION 


DATA 

NUMBER  GF  POINTS 

UCAH 

DIFFERENCE 

COEFFICIENT 

DIFFERCNCE 

CORRELATION 

COLUMN 

PC*  COtU**f4/AGWS 

(M) 

IN 

OF  VARIATION 

COEFFICIENT 

COMMENTS 

identification 
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'AM) 

(C’Ji 

OF  VARIATION 
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CODE 

TABLE  I 

iACV) 

SET  1 

wvhb 

30 

/IS-  20, 

m 

Wm 

49.1 

jmi 

0.9SC* 

woae 

WM8B 

22-23, ?», 

'  1 
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0.361* 

WP»B 

29,32-34,57-40, 

ypg  ^ 

LOW  f  ATTRIBUTABLE  TO 

WD8B 

42-  41,  AND 

33. i 

2.2 

0.934* 

FROM  1  TO  9  000  POINTS 

WP96 

umsb 

MDBS 

MM6G 

MP8B 

MO  SB 

52-6 

Si  .6 

30.3 

30.9 

30.3 

29.6 

30.  r 

m 

IN  EACH  COLUMN  PAIR 

OF  SO  POINTS. DISPERSION 

CONSISTENT  EXCEPT  FOR 

THESE  FEW  000  POINTS. 

(SEE  FIGURES  3  AND  4} 

0.6 

2.9 

0.990 

0.7 

2.2 

0.971* 

4.9 

0.9B9* 

mpbb 

29.6 

S3.0 

_ 

19 

IS 

1 

■9 

■■ 

0.999 

/JT-29,29, 

32-34, 

sb.i 

KB 

1 

0.994 

37-40,  AND  42 

3B.4 

39. • 

mm 

0.999 

i  T^^m 

36.4 

46.9 

27.1 

3.4* 

70.9 

9.6* 

0.971* 

■ISM 

30.9 

60.9 

DISTINCTION  DUE  TO 

1  VERY  ODD  POINT  OF 

MM6W 

27.1 

E  | 

70.9 

8.9* 

0.991 

18  POINTS  TOTAL 

MFbW 

29.1 

63.7 

MDSW 

30.5 

■ 

60.9 

MI'S* 

29.1 

mm 

S3. 7 

2.9 

0.972* 

TWO  VERY  ODD  POINTS 
APPARENT 

s 

13 

■  ■ 

/15-IS,  3?, 41-44, 

42.0 

44.3 

3.3 

0.349* 

W04B 

47,  4S,  AND  33 -3S 

48.6 

41  .0 

WM4B 

42.0 

M 

44.3 

0.9 

0.979* 

LOW  f  ATTRIBUTABLE  TO 

WP4B 

W04B 

39. S 

41.4 

41 .0 

2  TO  9  ODD  POINTS  IN 

ALL  THREE  PAINS. 

43. S 

WP4B 

39.6 

II 

43.4 

2.4 

0.929* 

SET  4 

IS 

MM4B 

/16,‘B, 26-29, 

39.2 

0.4 

39.9 

0.990* 

MD4B 

34,41-46,  AND 

39.6 

54.C 

\  ONE  OOD  POINT  IN 

-'S-SS 

(COLUMN  MM4B  WAS  THE 

MM4B 

39.2 

99.9 

l CAUSE  OF  LOW  f  IN  BOTH 

MP4B 

39.1 

0.1 

94J 
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<  39.6 

39.1 

0.9 

94.2 

94.6 

m 

0.999 

MP4B 
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AFAPL-TR-70-23 


With  uncorrected  data,  AM  was  significantly  large  with  1  of  the  18  pairs, 
and  ACV  was  significantly  large  with  3  of  18  pairs.  The  correlation  coefficient 
was  significantly  iow  with  12  of  18  pairs. 

The  first  set  contained  six  column  pairs, each  with  30  points  per  column. 
Five  of  the  six  pairs  displayed  significantly  low  r,  but  this  lack  of  correlation 
was  attributable  to  the  existence  of  1  to  5  odd  points  per  column  pair.  Figures 
3  and  4  show  two  of  the  correlation  plots  in  question.  Correction  of  the  odd 
points  resulted  in  all  correlation  coefficients  being  greater  than  0. 990. 

Where  low  correlation  appeared  in  the  other  three  sets,  it  was  also  attrib¬ 
utable  to  from  1  to  4  odd  points  in  each  column  pair.  Correction  resulted  in  r 
being  greater  than  0.  990  in  all  cases. 

There  are  no  trends  evident  in  the  Am  and  ACV  columns  of  Table  IV. 
Choice  among  the  three  refiectometers  did  not  appear  to  influence  the  resultant 
SN ;  no  one  of  the  three  displayed  superiority. 

6.  EFFECT  OF  REFLECTOMETER  BACKGROUND  SHADE 

The  results  of  this  series  with  10  column  pairs  arranged  into  three  sets 
are  shown  in  Table  V, 

The  effects  of  refleetometer  background  shade  and  filter  medium  choice  are 
closely  coupled.  It  was,  generally  evident  that  refleetometer  background  shade 
had  no  significant  influence  when  used  with  Millipore  filter  medium,  but  had 
significant  influence  when  used  with  Whatman. 

Of  the  five  Whatman  medium  pairs,  four  displayed  significantly  high  AM, 
all  five  displayed  significantly  high  ACV,  and  four  displayed  signifi'”"*’ 
r.  The  contrary  was  true  with  the  five  Millipore  pairs.  None  of  the  Am  ox 
A  CV  was  significantly  large.  One  of  the  pairs  displayed  r  less  than  0. 990,  but 
this  was  attributable  to  one  odd  point  out  of  17  per  column  in  the  pair. 

The  magnitude  of  the  effect  with  Whatman  paper  was  markedly  displayed 
in  the  plots.  Two  of  the  five  plots  are  included  here  as  Figures  5  and  6.  All 
five  Whatman  plots  displayed  translated  functions,  indicating  that  the  uncor¬ 
rected  data  correlation  coefficients  were  deceptively  high. 
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TABLE  3T -RESULTS  OF  REF LECTOMETER  BACKGROUND 
SHADE  INFLUENCE  INVESTIGATION 


Data  Column 

I  d«n  t  it  ication 

Cod* 

Number  of 
Points  p«r 
Column/Rows 
:from  Tablel 
Inclusive 

Mean 

(M) 

Difference 

in  Means 

(Am) 

Coefficient 

of 

Variation 

(  CV  ) 

Difference  in 

Coefficients 
of  Variation 

(Acv) 

Correlation 

Coeff  ic  i  en  t 

(r) 

Comments 

Set  1 

Ij  1  \ 

W  MSB 

35 
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0  9  82* 
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43  6 

■ 
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.. .  V  1 
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MM  88 
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33. 1 

\  4 

2  7 

0  99  1 
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MD  8W 

34  5 

it  corrected 

Set  2 

I 

WM4B 

15 

42  5 

7.4* 

38  5 

8  5* 

0  989* 

r  deceptive!  y  h  igh 

WM4W 

/  IS, 26-28, 

49  9 

30  0 

( 

Plots  showed  both 

34,  41  44, 

/ 

correlation  liras 

WD48 

42  5 

3  3 

37  9 

4,0 

0.  96  3  * 

i 

to  be  very  transloiad 

WD4W 

45  8 

33  9 

53-56 

MM4B 

39.  0 

0  3 

50  2 

0  4 

0  997 

MM4W 

39  3 

49.  8 

M04B 

38  1 

1  .  1 

49  2 

1  .6 

0  99  9 

NIL  4W 

3L  2 

47  6 

Set  3 

♦ 

♦ 

WP8B 

17 

29  2 

6  9 

46  5 

7  8 

0.  9  94 

r  deceptively  high 

WP  BW 

/ 18  -2  0 , 

36  1 

38  7 

Plot  showed  obvious 

22  25, 

Ira  r.slat  ion  (Figure  6) 

MP8B 

32  34  , 

26  0 

2  5 

64  8 

4  3 

0  9  86* 

One  odd  point 

37-40 , 

MP8W 

and  4  2 

27  5 

60.5 

- - 

r  >  0  990  if  corrected 

Significonl  distinction  in  the  uncorrtcted  data  S •  <  criteria  in  Section  HI,  Pjrograph  2 
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Figure  6.  Correlation  Plot  of  WP8B-WP8W  Showli^  Influence  of  Background 
Shade  with  Whitman  Medium 
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The  white  background  produced  higher  averages  (M)  and  lower  dispersion 
per  unit  mean  (CV)  with  Whatman  medium.  White  appeared  to  be  superior  when 
used  with  the  Whatman  medium. 

With  the  Millipore  medium,  background  did  not  have  significant  influence 
on  the  results,  so  there  is  no  superiority  of  one  background  value  over  the 
other. 

7.  EFFECT  OF  FILTER  MEDIUM  CHOICE 

Table  VI  gives  the  results  of  this  series  with  10  column  pairs  arranged 
into  four  sets. 

The  AM,  A  CV,  and  r  indicated  significant  difference  in  results  taken 
with  Millipore  vs  Whatman  media. 

The  coupling  of  reflectometer  background  shade  and  filter  medium  choice 
effects  is  also  very  noticeable  in  this  series.  With  white  background,  Whatman 
gave  significantly  higher  averages  of  6. 1  to  11.  0  with  all  five  pairs.  With  the 
five  black  background  column  pairs,  the  Whatman  mean  is  higher  than  the 
Mi” 'pore  mean  by  2.  5  to  4.  1,  but  these  A  M's  appeared  to  be  significantly 
large  in  only  two  of  the  five  cases. 

Not  only  were  the  differences  large  overall,  but  the  correlation  plots 
revealed  that  the  differences  in  results  from  Whatman  versus  Millipore  media 
were  consistently  greatest  in  the  important  region  of  C'N  15  to  45.  Figures  7 
and  8  show  two  of  these  correlation  plots. 

Whatman  displayed  consistent  superiority  over  Millipore  medium  in  ail 
10  comparisons.  The  Whatman  column  means  were  hlgb<.?t  in  ail  10  cases, 
and  the  amount  of  dispersion  per  unit  mean  (CV)  is  lowest  for  Whatman  in  all 
10  cases.  Ba<  ‘ground  choice  affects  the  magnitude  of  this  superiority.  Whatman 
was  much  more  superior  to  Millipore  on  the  whtte  background.  The  same  trend 
was  consistently  evident  with  black  background,  although  the  magnitude  of 
Whatman’s  superiority  was  less  than  that  of  Millipore. 
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S«»  crit«rio  in  Section  m,  Porjgroph 
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Figure  8. 
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Correlation  Plot  of  WM4W-MM4W  Showing  Greatest  Difference 
Between  Results  with  Two  Filter  Media  at  Low  SN 
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SECTION  IV 

DISCUSSION  OF  RESULTS 

1.  SAMPLING  FLOW  RATE  CONSIDERATIONS 

The  analysis  indicated  that  doubling  the  sampling  flow  rate  (from  0.  0041 
to  0.  0085  sefs)  produced  a  reduction  of  2%  to  9%  in  S.  .. 

In  1954,  Watson  reported  that  the  need  for  isokinetic  sampling  (i.  e. . 
matching  the  sampling  velocity  at  the  probe  entrance  to  the  surrounding  stream 
velocity)  became  greater  as  particle  size  increased  (Reference  4).  Recent  work 
has  indicated  that  the  particulate  matter  in  aero  engine  exhaust  is  of  such  small 
size  as  to  make  the  need  for  isokinetic  sampling  superfluous  (References  1,  2, 
3,  and  5).  The  results  of  this  analysis  tend  to  corroborate  that  recent  work. 

The  small  sampling  flow-rate  effect  must  be  considered,  but  ir  does  not  appear 
to  be  large  enough  to  justify  the  complexity  and  effort  involved  in  isokinetic 
sampling.  Merely  specifying  a  standard  flow-rate  value  seems  to  be  proper 
and  sufficient. 

Some  smoke  measuring  systems  employ  sample  volume  and  sampling  time 
measurements  to  determine  ilow  rate.  The  analysis  also  indicated  that  such 
a  more  precise  yet  laborious  procedure  for  determining  flow  rate  is  super¬ 
fluous,  The  analysis  tends  to  indicate  that  variations  in  flow  rate  of  as  much 
as  10$  will  produce  variation  in  results  (SN)  of  less  than  1%. 

2.  CHOICE  OF  REFLECTOMETEH 

The  analysis  indicated  that  all  three  reflectomoters  produced  substantially 
the  same  results.  The  differences  in  data  column  means  were  2$,  to  11%,  but 
there  was  no  evidence  of  any  one  reflectometer  producing  superior  quality 
results. 

3.  CHOICE  OF  FILTER  MEDIUM  AND  REFLECTOMETER  BACKGROUND 
SHADE 

The  results  showed  that  fjR  -  media  and  background  shade  effects  were 
coupled. 
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Combinations  of  filter  media  and  reflectometer  background  shades  are 
ranked  in  '’'able  VII.  The  M  and  CV  averages  in  Table  VII  were  prepared  from 
T  ables  V  and  VI.  The  corresponding  SD  averages  were  prepared  from  computer 
calculated  SD  not  reproduced  in  this  report. 

TABLE  SHI  SUPERIORITY  RANKING  OF  FILTER 


MEDIUM -REFLECTOMETER  BACKGROUND  SHADE  COMBINATIONS 


Comb  in  otion 

Average  Mean 

(Maximize  for 

Superiority ) 

Average  CV 

(Minimize  for 

Superiority ! 

Average 

Standard 

D  tv  ia  tion 

Ranking 

Whotmon/whit* 

43.  0 

37.  i 

15.4 

Best  Combination 

Who tmon/blnck 

36.  9 

44.  6 

16.  1 

Mi  1  i  i  pc-  r3  /  «*  h  i  t  • 

34.3 

56.4 

18.3 

)  Either  Combination 

M  I II 1  port  /black 

33.  4 

55.9 

18.1 

[  Least  Desirable 

'  (insignificant 

Difference  Between 
These  Two) 

Whatman  with  white  background  displays  significant  superiority  (highest  M, 
lowest  CV)  over  all  other  combinations. 

With  Millipore,  the  differences  between  results  on  either  white  or  black 
background  are  not  significant. 

The  combined  results  given  in  Table  VII  are  for  the  full  range  of  smoke 
levels  investigated  during  the  tests  (approximately  SN  of  5  to  70).  However,  the 
correlation  plots  revealed  that  the  distinction  between  Whatman  and  Millipore 
media  is  even  greater  in  the  most  important  region  of  SN  from  about  10  to  45 
(see  Figures  7  and  8). 

The  effect  of  black  background  is  to  decrease  the  magnitude  of  Whatman 
superiority  by  decreasing  the  overall  average  SN  from  43.  0  on  white  to  36.  9 
on  black  (Table  VII).  It  is  significant  to  note  that  the  higher  CV  of  Whatman  on 
black  versus  Whatman  on  white  is  largely  attributable  to  this  reduction  in  SN ; 
the  amount  of  dispersion  in  Whatman  data  is  about  the  same  with  both  black  and 
white  background  shades.  Conversely,  the  even  higher  CV  averages  of  both 
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Millipore  combinations  are  primarily  attributable  to  more  dispersion  (higher 
standard  deviations),  in  addition  to  lower  column  means. 

The  superiority  of  Whatman  paper  has  been  previously  implied  if  not 
explicitly  denoted.  Bagnetto  (Reference  2)  evaluated  three  smoke  measuring 
systems  and  concluded  that  the  Von  Brand  system,  which  uses  Whatman  No.  4 
medium,  was  significantly  superior  to  the  AED  system  that  used  Millipore. 
(The  third  system,  the  B.  P.  Hartridge  nonfiltration  type  based  on  light 
absorption,  was  ranked  slightly  above  the  AED  system,  yet  still  significantly 
below  the  filtration  type  system  using  Whatman  medium. ) 

It  may  be  possible  to  reconcile  the  differences  in  results  obtained  on 
Millipore  versus  Whatman  media.  One  theoretically  possible  tack  for  making 
the  results  of  Millipore  medium  approximately  equal  to  those  of  Whatman 
medium  is  discussed  in  Appendix  V. 
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SECTION  V 

CONCLUSIONS  AND  RECOMMENDATIONS 

1.  CONCLUSIONS 

The  subject  data  can  be  used  to  draw  statistically  valid  conclusions  about 
four  parameters:  sampling  flow  rate,  choice  of  reflectometer,  choice  of  filter 
medium,  and  reflectometer  background  shade. 

Sampling  flow  rate  had  a  small,  yet  consistent  and  significant,  influence. 
The  higher  flow  rate  tested  (0. 0085  scfs)  produced  SN  results  that  were  2%  to 
9%  lower  than  results  with  the  lower  flow  rate  (0. 0041  scfs).  Since  there  was 
not  enough  data  to  compare  results  at  each  of  the  four  engine  power  levels  used 
during  testing,  no  firm  statement  can  be  made  about  the  need  for  isokinetic 
sampling.  However,  the  analysis  does  tend  to  corroborate  previous  work  that 
concluded  that  the  isokinetic  sampling  requirement  is  superfluous  when  sampling 
exhaust  smoke  from  aircraft  gas  turbine  engines. 

All  three  retlectometers  used  to  rate  spots  (MacBeth,  Densichron,  and 
Photovolt)  were  found  to  produce  essentially  the  same  results.  No  one  of  the 
three  demonstrated  superiority. 

The  effects  of  filter  media  choice  and  reflectometer  background  shade  are 
closely  coupled.  Whatman  filter  medium  evaluated  on  the  white  background 
produced  the  best  results.  Whatman  filter  medium  on  black  background  gave 
significantly  lower  SN,  but  the  dispersion  of  data  with  this  combination  was  not 
significantly  different  from  Whatman  on  white.  Whatman  on  black  was  the 
second  best  combination. 

Background  shade  did  not  significantly  influence  results  obtained  with 
Millipore  filter  medium.  Results  of  Millipore  with  either  background  were 
significantly  more  dispersed  than  the  results  of  Whatman  with  either  back¬ 
ground. 

Whatman  medium  with  either  background  was  superior  to  Millipore  with 
either  background. 
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2.  RECOMMENDATIONS 

The  SAE  system  for  measuring  aero  engine  exhaust  smoke  should  specify 
a  single  sampling  flow-rate  value.  It  should  also  specify  a  rotameter  or  other 
simple  device  for  direct  measurement  of  flow  rate.  The  influence  of  sampling 
flow  rate  on  results  is  not  great  enough  to  justify  the  need  for  isokinetic 
sampling,  nor  is  the  influence  great  enough  to  justify  more  precise,  indirect, 
means  of  determining  sampling  flow  rate. 

The  SAE  document  should  either  specify  use  of  any  of  the  three  reflectom- 
eters  tested,  or  otherwise  ensure  that  an  inferior  instrument  is  not  allowed. 

The  use  of  Whatman  No.  4  filter  paper  and  black  reflectometer  background 
(i.  e. ,  absolute  background  reflectance  of  5%  or  less)  should  be  specified.  The 
analysis  determined  that  white  background  was  superior,  but  its  use  is  un¬ 
desirable.  There  is  some  evidence  that  white  background  tends  to  exaggerate 
the  differences  ir '  veen  values  of  all  parameters,  as  the  analysis  indicates  it 
does  with  filter  media.  In  the  first  set  of  Table  VI,  the  ACV  of  columns 
WM8B  and  WD8B  is  virtually  identical:  47.  7  -  47,  6  =  0. 1.  Yet  the  ACV  of 
WM8W  and  WD8W  is  much  larger:  40.  3  -  36. 8  =  3.  5.  Since  it  is  impossible  to 
specify  limitations  on  aii  parameters,  the  use  of  white  background  could  cause 
poorer  results  by  affecting  unspecified  parameters.  The  analysis  indicates  that 
use  of  black  background  does  not  appreciably  increase  data  dispersion,  although 
its  use  will  cause  a  slight  loss  in  precision  due  to  the  absolute  value  of  SN  being 
lower.  This,  in  the  author's  opinion,  is  a  justifiable  tradeoff  to  guard  against 
unknown  factors. 
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APPENDIX  I 

SUMMARY  OF  TEST  PROGRAM  AND  MEASUREMENT  SYSTEM  EVALUATED 

By  June  of  1969,  SAE  Committee  E-31  kad  decided  that  the  standard  mea¬ 
surement  system  should  be  an  indirect,  filtration  type  not  unlike  most  systems 
in  use  today.  A  test  program  was  conducted  by  the  Committee  at  the  Federal 
Aviation  Administration's  experimental  center  in  Atlantic  City,  New  Jersey. 

The  program  was  intended  to  experimentally  examine  the  tentative  smoke  mea¬ 
surement  system. 

A  J-57  turbojet  engine  was  used  to  generate  smoke.  The  elements  of  the 
measurement  system  being  evaluated  were  provided  by  various  Committee 
members. 

The  design  of  the  experimental  program  is  unknown  to  this  author,  although 
it  is  known  that  the  testing  sequence  used  did  pproach  being  random.  Data  w  s 
obtained  with  different  combinations  of  values  of  the  following: 

•  engine  power  level  (4  values) 

•  filter  media  (2  types) 

•  filter  media  holder  (2  types) 

«  sample  size  {standard  volume  of  exhaust  gas;  4  valuer  for  each 
of  the  2  types  of  filter  media) 

•  sampling  flow  rate  (2  values) 

•  sampling  probe  angular  orientation  with  respect  to  the  direction 
of  engine  exhaust  gas  flow  (3  values) 

•  sampling  probe  position  along  the  engine  exhaust  gas  b.  ream 
flow  path  (2  values) 

•  sampling  line  length  (2  values) 

•  sampling  line  size  (diameter  -  2  values) 

•  rumpling  line  material  (2  types) 

•  sampling  line  temperature  (2  values) 

Data  reduction  introduced  two  more  parameters  with  variable  values: 

•  reflectometer  (3  types) 

•  reflecionMer  background  shade  (2  values) 
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Not  all  parameter  values  were  changed  for  all  runs.  Yet  there  was  enough 
change  so  that  no  two  rows  of  Table  I  contained  all  the  same  parameter  values. 

Mor°  than  200  data  points  were  taken. 

The  basic  configuration  of  the  measurement  system  and  operating  pro¬ 
cedure  were  the  same  throughout  testing.  A  given  sample  size  was  drawn  at 
a  given  flow  rate  from  the  engine  exhaust  through  the  sampling  probe,  sampling 
line,  and  filter  media  holder  with  a  vacuum  pump.  A  rotameter,  positive  dis¬ 
placement  volume  space  meter,  and  pressure  and  temperature  gauges  were  used 
downstream  of  the  pump  to  measure  the  sample  before  it  was  discharged  to 
the  atmosphere.  The  flow  time  of  each  sample  was  also  measured.  The  system 
was  heated  throughout  testing.  A  filter  holder  bypass  line  was  used  to  maintain 
flow  rate  in  the  system  when  a  sample  was  not  being  taken. 
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APPENDIX  II 

DATA  REDUCTION  PROCEDURES 


1.  SAE  GROUP 

The  smoke  spots  were  read  with  meter-background  combinations  to  rate 
them  in  terms  of  either  absolute  reflectance  or  optical  density,  depending  on 
which  meter  was  being  used.  These  readings  were  then  used  to  calculate  SN. 
The  definition  of  SN,  the  dimensionless  term  used  to  quantify  smoke  emission, 
is: 

SN  =  10 0  \l  -  —  '  ,  wh«r«  { 

Rs  -  absolute  reflectance  of  the  sample  spot 
Rw  -  absolute  reflectance  of  clean  filter  media 

The  relationship  between  optical  density  (OD)  and  absolute  reflectance  (R) 
is: 

=0  •  (^ ) 

The  SAE  data  reduction  group  used  graphs  combining  Equations  1  and  2  to 
obtain  SN  for  spots  used  in  terms  of  optical  density  (the  MacBeth  meter). 
Equation  1  was  used  to  calculate  SN  for  spots  read  in  terms  of  absolute  re¬ 
flectance  (the  Photovolt  meter). 

SN  will  vary  with  the  sample  size.  It  has  long  been  accepted  to  report  SN 
and  other  quantifiers  of  smoke  for  a  certain  sample  size  (standard  cubic  feet) 
per  unit  filter  medium  area  (square  inch).  This  quantity  is  termed  "Q.  " 

It  has  also  been  accepted  practice  to  use  a  specific  Q  value  dependent  on 
filter  medium  choice.  The  ''standard"  Q  for  Whatman  medium  is  0.  300  scf/sq 
in,  and  the  value  used  for  Millipore  is  0.  0565  scf/eq  in.  These  are  widely 
used, although  apparently  arbitrary,  values. 

The  SAE  Group  calculated  Q  values,  and  then  plotted  these  as  abscissa 
versus  the  corresponding  SN  as  ordinate  on  log-log  paper.  A  curve  waa  then 
fitted  to  these  points,  and  then  the  SN  values  were  read-off  for  Q  values  of 
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0.300  and  0.  0565  scf/sq  in  for  What  e  andMillipore  filter  media, 
respectively.  These  SN  are  the  values  reported  in  Table  I. 

2.  GE  GROUP 

The  GE  Group's  procedure  differed  from  that  of  the  SAE  Group  in  the 
manner  in  which  the  effect  of  different  sample  sizes  was  weighed.  The  GE  team 
has  a  standard  smoke-spot  data-reduction  routine  based  on  the  use  of  "loading 
curves"  for  Whatman  and  Millipore  media.  The  loading  curve  is  a  plot  of  micro¬ 
grams  of  carbon  (smoke  particulate  matter)  as  abscissa  versus  optical  density 
of  the  resultant  spot  as  ordinate  (Figure  9).  Spot  size  is  a  necessary  parameter 
of  such  curves. 

The  GE  data  reduction  routine  was  completely  computerized.  All  raw  data 
to  compute  SN  and  Q  were  input  (reflectometer  readings,  sample  size,  etc. ). 

The  routine  calculated  SN  and  corresponding  Q,  and  then  based  on  leading 
curve  factors,  "corrected"  each  SN  to  the  "proper"  value  for  Q  of  0.  300  or 
0.  0565  scf/sq  in,  depending  on  which  filter  media  was  used.  The  resultant  SN's 
from  this  procedure  are  those  listed  in  Table  I. 
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APPENDIX  III 

ANALYSIS  COMPUTER  PROGRAM  Cr  .wR.El 

The  CORRE1  calculation  routine  is  written  in  FORTRAN  IV,  Version  13 
(Reference  10).  The  short  main  routine  MP"  calls  the  primary  subroutine 
"MAIN.  "  MAIN  in  turn  calls  the  subroutines  "CORRE"  and  "GP.  "  The  final 
subroutine  "DATA"  is  a  short  dummy  element  used  only  to  satisfy  a  call  from 
CORRE.  (DATA  is  included  to  avoid  having  to  modify  CORRE. )  Subroutines 
CORRE  andGP  were  taken  from  References  7  and  8,  icspectively. 
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MAIN  ROUTINE  MP 


(.•VALUATION  ur-  COMM.  E-31  TEST  DATA  -  ROUTINE  CURREi 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

r 

L* 

c 

c 

r 

c 

c 

c 

c 


INPUT  -  MT  -  NUMBER  or-  COL  ONS  UE  l)aTu  (12  MAX. 

PER  CASE) . 

-  NT  -  ROWS  OF  DATA  PER  COLUMN.  b'}  MAX  .  AN!) 

MUST  NOT  BE  LESS  THAN  r-iT . 

-  NC  —  NUMBER  OF  CASES.  EACH  SET  OF  12  OR  EEi'ER 

COLUMNS  t(JK  EACH  DIFFERENT  i\*T  VALIjc  t 
CONSTITUTE  A  CASE. 

-  JP  -  NUMBER  UE  CULUMV  PAIRS  TO  BE  C: :  R  R  E  L  A  T  Eu . 

-  SNPCU!  -  ALPHAMERIC  COLUMN  COPE  (6  SPACES  --AX.  ). 

J  IS  THE  COLUMN  »NUEX  NUMBER. 

-  JX (  I  J)  »  J?  (  I  J) 

-  INDICES  UP  CULUMI-  ’•'$  TO  dE  CORRELATE!;. 

I J  =  1  IS  THE  FIRST  PaIR  ,  ETC . 

-  SN(ItJ)  -  DATA,  INPUT  COLUMN  BY  COLOR*'. 

-  H I  (  I  )  —  A  42  '■rACE  MESSAGE  UE  ivHICH  r  >  /v  y  4  RCivS 

WE'  \  PUT .  I  IS  A  DUmhY  SUBSCRIPT. 
OUTPUT  -  DATA,  COLUMN  H.  JLUmiv,  WITH  CULUmiv  IDENTIFY  IF  G 
CODE  AnD  STATEMENT  OF  r-hlCH  ROWS  RERe  I  "'PUT . 

-  MEANS, so, STANDARDIZED  MEANS  AND  SO,  «ND  C(  iE EE  1 C 1 E NT S 
OF  VAfi  I  AT  I  (Tim  ,  EUR  ALL  INPUT  COL  On  n  S . 

-  CORRELATION  COEFFICIENTS  FOR  SELECTED  COLD'-- 
PAIRS. 

-  PLOTTING  OE  EACH  CORRELATED 
COLUMN  PAIR. 

i )  I  m E n S I  UN  Y(lDO) 

RE  AD (b, 10)  NC 
IC  =  1 

BO  K E  AD ( b ,  10  )  NT 
CALL  MAIN  ( Y , N  T ) 

IC  *  IC+I 

I E ( I C-nC )  bO , bO , 60 
10  EURMAT(Ib) 

60  STOP 
End 
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PRIMARY  SUBROUTINE  MAIN 


Si  Ih*  OUT  I  >■ ■  r  (  Y  ,  n  T  ) 

o I  SI  I11':  XhAK  (  12)  ,  STO(  ]  °)  ,kX(  144  i  ,RI78)  ,tJ{  12i  ,0<  12  )  ,T  i  12!  ,AMNi  12  J 
i,S'i(i£|,  SAi-  (  12  )  ,  3SO(  12,  ,  JX(  20  )  ,  JZi  20  I  ,  SNPC  ( 1 2  >  *  C  V  (  1  i>  ’ 

2 S (SO  ,  12  !  .X  (600  )  >Y  '.NT,  2  ;  ,A  (2) ,  PLOT!  36,60)  ,KP  (20  )  .KI  (7) 

kfcAi)(  b  ,  10  )  N'T,  JP  ,  (SNPCI  J),J=1,HT) 

id  (  jy.  ( I  j) ,  j?  (1  j  ) ,  I  j  =  l ,  j p ) 

-Nh  Ai'O  ,  i2  )  (  (  SN  !  I  ,  J  )  ,1  =1,M1 ,  J=l,  HI ) 

KP  AO  (  b,  22)  (PI(!),  1  =  1,7) 

COMPUTE  t-  «('*  S ,  STANDARD  D  tV  I A  T  IUnS ,  A  NO  CO  fc  EE  IC  I  ENTS  UP  VARIATION 


J  i>  =  =  i 

i'll  AO  J=l,,« 

’  So  i  =  1  ,1- 
J  S  =  J  S+  1 

31'  *{JS)  =  Sw(I  ,J) 

A i 1  Ciil>  T  iNDt 

Call  COKkfc  (N  ,n,lU,X,XBAR»STD,kX,k,h,U»T) 

Oil  70  J  =  1 , 

C V ( J !  =  lOO.vSTHi J ) /XbAk ( J  ) 

A'-'v(J)  =  XhAK(J) 

7o  SD(J)  =  STo(J) 

CUi-PUTE  AVbkAGfc  MEANS  AND  AVFBAGfc  STANDARD  DEVIATIONS,  AND  OUTPUT 
wITh  I '  POT  DATA  AND  MEANS,  STANDARD  DEVIATIONS  AND  COEEP.  UP  VARIATION 

iv  =  1*1  T 

'•i  =  > 

I  S=o 

'  H  i  HO  I  =  1  ,  N 
I S= I  5+1 

P ( 1  X ( I  S)  =  XhAK  (  I  ) 
no  ns  I  =  1,n 
I S= I $+  1 

PS  X ( I S )  =ST»( I  ) 

1,  a L  !..  ^  0  Kk.  t  (  o  2 U  ,  X  ,  X  BAR  ,  STD,kX,  K  ,  R  ,  I) ,  T  ) 

AAi-iN  =  X  BAK  (  1  ) 

A  Si )  =  X(>AK  ( l  ) 
i*.  =  h  T 

Dll  40  J  =  1  ,  1*1  T 
SAi'i(J)  -•  Ar-ilv  (  J  )  /AAl'IN 
40  SSo(J)  =  Sl>(J)/ASD 

-H I Tt (6, 14 )  ( SNPC( J) , J= 1 , NT ) 

DO  ob  1=1, NT 

4b  OKI  TP  (6,  IP  ,  (  Sim  ( 1  ,  J  )  ,  .)=  1 ,  oT ) 

wk I  TP (h,23)  ( k  2  ' I  ),  1  =  1,7) 

WK  I  T  H  ( t> ,  1 3 ) 

|v  =  IV  T 

Dll  100  J  =  1 ,  i*i T 

inn  WK  1  11(0,14)  SNPC  (  J  !  ,  N ,  AMIm (  J  )  ,  SD(  J  )  ,  SAM  (  J  )  ,  SSI)  (  J  )  , C  V  (  J  ) 

n  K  I  T  B  (  6 , 1  *  '  A  A  i'i  N ,  A  S I ) 
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COMPUTE  AND  OUTPUT  CORRELATION  CDEFEIC  IEnTS 

WR I T  E ( 6 ,  1  b ) 

=  2 
i  J=  1 

no  is  =  o 

JI  =  JX(IJ) 

J2  = JZ ( I J ) 

DO  12U  1=1,  iv. 

IS=IS+1 

120  X  IIS!  =  Sim  (I  ,Jli 
IH!  i25  I  =  l,i'' 

I S=I S+l 

12 5  XI  IS)  =  Soil  ,J2) 

CALL  CUKKE  (N,l-:f  IO,X,XBAk,STf>,kX,R,B,l)(T  ) 

NR1  Tfc  (6, 17  )  SNPCI  Jl  i  , SNPC!  J2  )  ,k  i  2  ) 

R  P  (  I  J  !  =  R  (  2  ) 

I  J=I J+l 

IF(IJ-Jp)  110,110,130 

X-Y  PLUTT I  iv  ft  OF  CORRELATED  COLUMNS 

130  C  OImT  I  nOE 

!  NT  blit  K  S  ,  w 
DATA  A/lH»,iH+/ 

N=2 

(.1  =  M  T 
L  =6 
S  =  30 
W  =  60 
Ln  =  36 
I  J=  1 

ISO  l-i  K  I  T  F  ((>,21) 

Jl  =  JX I  I J  ) 

J  2  =  j7  (  I  J  ) 

00  160  1  =  1,1' 

XU)  =  SN(I  ,  Jl) 
v  1 1  ,1)  =  Si"  1 1  ,  J  1 ) 

160  Y  (I  ,2)  =  SI"  (  I  ,  J2) 

CALL  G  P  (  X  ,  Y  ,  L  ,  S  ,  I" ,  m  ,  W ,  L  N ,  A  ,  P I.  O  T  ) 

MR  I TF (6,20)  SNPC(Jl)  ,  SNPCIJ2)  ,KP(  I  J  ) 

IJ  =  I J+l 

IF(IJ-JP)  150,150,200 
C 

c 

10  FORMAT ( 2  I  6 , /H ( 4X , A6 )  ) 

11  FORMAT (4(13,15, lOX ) ) 

12  FORMAT ( i6P5. 1 ) 

13  FORMAT! 1H1,30X,60H  AVtR AGES ,  STANDARD  DEVIATIONS  AND  COEFFICIENTS 
10E  V  AH  I  AT  I  Ol'  ///  1  5X  ,  7H  COLUMN, 

2  7X,7H  NOi-'BEK  ,  3X  ,  lit'  ARITHMETIC,  SX.SH  S  TANDAKD  ,5  X  ,  1 3H  STANDARDIZE 

3D , 6 X , 1 3H  STANDARDIZED, 4X,16h  COEFFICIENTS  0F/10X, 

A  1 6H  (PARAMETER  SET),2X,«H  OF  DATA,6X,5H  NFAN  ,  P.< ,  10H  DEVIATION, 

5  5  X  ,  1 1 H  ARITH  MEAN, 5X, 15H  STAN  DEV  I  AT  ION , AX , 15H  VARIATION  (CVI/26X 
6 ,  1 2H  !SN)  POINTS, 3X,6H  ( A  MM } ,  1 0  X , 5  H  (SD),hX,llH  ( AMN / A ANN  )  , 7 X , 9H  ( 
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7SD/AS0  )9X  ,  15H  (100  X  50/AMN)//) 

14  FORMAT! 16X,A6,8X,I2,5F16.3) 

15  FORMAT! 1H0,40X,34H  AVERAGE  ARITHMETIC  MEAN-  (  AAMN  )  = ,  F  10 . 3/40  X  ,  36  H 

1  AVERAGE  STANDARD  DEVIATION  (  A S 0 )  =,110.3//////) 

16  FORMAT (4BX.2H  X,10X,24H  CORRELATION  COE FF IC I  ENT/45  X , 2H  Y,14X,16H  b 
1ETWEEN  X  AND  Y//J 

17  FORMAT (42X ,A6/42X ,A6, 10X,F 10 .3,/ / ) 

18  FORMAT!  10X.12F10. 2) 

19  FORMAT  ! 1H 1 , 60X , 1 1H  INPUT  DATA// 10X , 1 2 [ 4X  , 46 ) / / ) 

20  FORMAT  (  1H0, BOX,  18H  ABSCISSA  (X)  IS  ,  46/5  1 X  ,  IfiH  UR!)  I NA  i  E  !Y)  IS 
1,A6/40X  543h  CORRELATION  COEFFICIENT  BETWEEN  X  AND  Y  = , F  10  -  5 ) 

21  FORMAT  (  1H  l ,  40X  ,  35H  PLOTTING  OF  CORRELATE!)  COLUMN  P  A  I  R // 30  X  ,  44H  PE 
1RFLCT  CORRELATION  LINE  IS  DRAWN  WITH  00 T S/30 X , 40H  ACTUAL  COLUMNS  P 

2  L  0  f  T  !  G  IS  WITH  CROSSES) 

22  Forma  f  !  7  a  i ) 

23  FORMAT  (  :’••(!,  JXth'-K  THIS  OATA  IS  FROM  THE  FOLLOWING  KUwS ,  INCLUSI VE 
i, READING  FROG  COLUMN  TOP  TO  BOTTOM  -  ,7A6) 

0 k  HTu R i'i 
END 
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SUBROUTINE  CORRE 


SUBROUTINE  CUKKE 
PlJkPUSt 

COMPUTE  MEANS, 


STANDARD  DEVIATIONS,  SUMS  UP  CROSS-PRODUCTS 


UP  DfcVlATIUMS,  AND  CORRELATION  CHEER IC IENTS . 
USAGfc 

Call  CUkkF  ( N, M , 1 0 , X, XBAK , STD , kX ,K ,b ,0 , T ) 

DESCk  I  PT!  UN  (IP  PARAMETERS 

IV  -  NUMBER  LIF  OBSERVATIONS, 

.-I  -  NUMBER  UP  VARIABLES. 

ID  -  OPTION  CODE  F0»  INPUT  DATA 


COkRPOO 1 
.  CL)  Rk  LOO  2 
CORK  POO  3 
CORK tOO 4 
CURRFOOS 
COKRPOO 6 
CORK  POO  7 
CORK  POOH 
Cl  I RR  POO  4 
CUKKEO  10 
CORK  P0  11 
CORREO  1  2 
COkR  EO 1 3 
CORK EO 14 
CORK  FO  1  5 
COR  RFC)  In 


0  IP  DATA  AKr  TO  BE  KFAU  IN  FROM  INPUT  DEVICE  IN  THFCUKR EO 1 7 


SPECIAL  SUBROUTINE  NAMED  DATA.  (SPE  SObROOTIMPS  CORREOIP 

USED  BY  THIS  SUBROUT  livp  bELOW.)  CORKED  19 

1  IP  ALL  DATA  ARP  ALREADY  1  iv  CORF,  CURKE020 

-  IP  10=0,  THE  VALUE  OF  X  IS  0,0.  C0RREU21 

IP  10=1,  X  IS  THE  INPUT  MATRIX  ( N  BY  M)  CUNT  A  IN  I  iv  G  C0kRE022 

DATA.  CURKP023 

-  OUTPUT  VECTUk  OP  LENGTH  h  CONTAINING  M  fc  A  in  S  •  C0KKE024 

-  UOTPOT  VECTUk  OP  LENGTH  h  CONTAINING  S T AND Ak I )  CURKF0  2  5 

L»F  V I  ATI  PINS  .  C0KKF026 

~  OUTPUT  MATRIX  (m  X  m)  COnTaIN.  >)G  SONS  OP  CROSS-  C0kkE027 

PRODUCTS  OP  DEVIATIONS  FROm  mFANS.  COkkP02P 

UOTPOT  MATRIX  (ONLY  UPPER  TRIANGULAR  PORTION  OP  THE  CURRF029 


SYMMETRIC  MATRIX  UP  M  BY  m)  CONTAINING  CukkPL AT  ION 
CnpPFICIENTS,  (STORAGE  NUDE  OF  1) 
rt  -  OUTPUT  VECTUk  OF  LENGTH  M  CONTAINING  TFifc  U I  AGON  a 

UP  THE  MATRIX  OP  SUNS  IIP  C  K(i  SS-  PRODUCT  S  UP 
DEVIATIONS  PROM  MtANS. 

U  -  WORKING  VECTOR  UP  LFNGTh  m. 

T  -  WORKING  VECTOR  OE  LENGTH  M. 

R  Pi-  AK  K  S 

n  MUST  BP  GREATER  THAN  OR  EOlJAL  TO  m. 

SUBROUTINE'S  AND  EUnCTIUN  SUBPROGRAMS  kFL'UIKFII 

MATA  (  M ,  I  > )  -  THIS  SUBRl'l  IT  I  l\'E  MUST  BE  PKUVIllEI)  BY  ThE  USER. 

(1)  IP  10=0,  THIS  SUBROUTINE  IS  EXPECTED  To 
FURNISH  aN  OBSERVATION  IN  VECTOR  I)  EkOM  A N 
EXTERNAL  INPUT  DEVICE. 

(2)  IE  10=1,  THIS  SUbRUUT  I NF  IS  Mill  USED  BY 
CUKKF  BUT  MUST  EXIST  In  JOB  DECK.  IE  IlSPk 
HAS  MOT  SUPPLIED  A  SUBROUTINE  named  DATa, 
THE  FOLLOWING  IS  SuGGPSTFD. 

SUBRUUT I  Nr  DATA 
RETURN 
E  m  1 1 

I-1  E  T HI  III 

PKHUUCT-hUmEnT  CORRELATION  COEFFICIENTS  ARE  Cni'ipuTEM. 


CORKED  30 
CO  PR  ED  31 
ClI  k  K  E  0  32 
ClJP.K  R)  33 
CORK EO 34 
CORK  1-0  3S 
Cl  i  R  R  El)  3  H 
Cork FO 37 
CIIRKE0  38 
I'.i  i  PR  Ell  3 vi 
LI  I  P  R  EO  4|) 
CORREO 41 
CORKE042 
Cork  RO  4'j 
(.UK  PRO  44 
Cl  IkK  El)  4S 

Cl  I R  R  E  0  4  E 
Cl  IRKED  4  7 
I'.I  l  p  K  R0  4H 
Co  KR  El  I  44 
CORK  Eli  So 
LURK  Ell  SI 
CORK  FOB? 
l,l  IRK  Fl)  s  ; 

Cl  I  KK  Ell  S4 
Cl  'RK  Rll  y< 
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c 

c 

c 


c 

c 


c 

c 

c 

c 


c 


SimkUuTIiNfc  Cl  IKK  e  (N,M,IO,X,XBAR,STO,RX,K,B,D,T) 
DImFivSIdN  X  (  1)  ,XBAR  (  1>  ,STD(  1  )  ,kX(  1)  ,R(  1)  ,h(l)  ,D(  1)  ,T<  1  ) 


IF  A  IHHJBLt  PRECISION  VERSION  OF  THIS  ROUTINE  IS  DESIRED,  THE 
C  I.'i  COLUMN  I  SHOULD  Bt  RFHJVEI)  FROM  THE  DOUBLE  PRECISION 
STATEMENT  WHICH  FOLLOWS. 

DOUBLE  PRECISION  X B AR , S Tl 1 , K X  ,  R  ,  B  ,  T 

ThE  C  i"llST  ALSU  BF  REMJVFH  FRO M  DOUBLE  PRECISION  STATEMENTS 
APPEARING  In  UTHFR  RLHJTINFS  USED  IN  CONJUNCTION  WITH  THIS 
RUUTI I'.'E  . 


THE  DOUBLE  PRECISION  VERSION  OF  THIS  SUBROUTINE  MUST  ALSO 
COmTaIn  DOUBLE  PRECISION  FORTRAN  FO. „TIOWS.  SORT  AND  ABS  In 
STATEMENT  220  MUST  BE  CHANGED  TU  DSORT  AND  DAKS. 


INITIALIZATION 

DO  100  J=l,h 
B  (  J ) =0  .0 
1(H)  T(J)=D.O 

K-  Im«-  ■+n  )  /2 
011  10  2  I  =  1,K 
102  R ( I  ) =0  .0 
Em  =  I" 

L  =0 

I F ( 1 1 1 )  ins,  127,  105 

DATa  ARE  ALkbADY  IN  CORE 

105  DO  108  J=1,N 
00  10  7  I  -  1 , N 
L  =L+  1 

107  T(J)=T(J)+X(L) 

X  B  AR  (  j  )  =  T  (  J  ) 

108  T ( J ) = 1 ( J ) /En 

00  115  1  =  1  ,  N 
J  K  =  0 
L  =  1  -i'1 

00  1 10  J=1,M 
I  -L  +i' 

0(J)=X(L)-T(J) 

110  B(J)=r(J)+0(J) 

I'll  in  J=  1  ,  I-I 

oi '  in  k=i,j 

J  K  s  J  K  +  I 


CO  R  R  FO  5  6 
CORKED  5  7 
COR RED  5 8 
CLIRREO  59 
CO  RR  FO  60 
CORK  EO  61 
CORK  FOB  2 
CORKED  63 
COR RED  64 
CORK  EO  65 
CORK FO 66 
CO  RR  EO  67 
C0RKF068 
CORK EO 69 
CORKED  70 
CORK EO 71 
Cl )  RR  EO  72 
CORK FO 73 
CO  R  R  EO  74 
CORKED  7 5 
CUR RED  7 6 
CORK  EO  77 
CLIRREO  7  8 
CORK  FO  79 
CORK F080 
CO  kK  rO  8  1 
CORK  FOE  2 
CORK  EO  8  3 
CORK FO  84 
CORK  E08  5 
CORKED 8 6 
CORK  EO  87 
CORK FO  88 
CORK  F089 
CORK  F090 
CORK EO 91 
CORKED 9 2 
CORK EO  9  3 
CORKED 9 4 
CORK  EO 9 5 
C0RRE096 
CORKED  9  7 
CO  RR  FO  98 
CURK  EO  99 
CURRF  100 
CORK  F  10  1 
Cork  e  it)  ? 
CORK  E  JO  3 
CORK E  10 4 
CORK  E  10  5 
Cf  I  kK  E  10  b 
CORK  E 10 7 
CORK E  10  8 
CO  KK  E  10  9 
CORK!  1  10 
CORK  Fill 
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11D  k( JK)=K( JK)+U( J)#D(K( 

GO  TO  20b 
C 

C  K  6  AO  UBSbkVAT!  OMS  AND  CALCULATE  T Em RU  K  AR  Y 

C  I'lfcANS  FROn  THbSE  DATA  IN  T(J) 

C 

127  IF(N-h)  130,  130,  135 
130  KK.=D 

GO  TO  137 
133  KK=M 

137  l)U  140  1  =  1, KK 
CALL  I)aTA  (n,0) 

Du  140  J  =  1 ,  i" 

T  (  J  )  =  T ! J ) +0 ( J ) 

L  =  L  +  1 

140  XX (L  )=D( J) 

FKK  =KK 

DO  130  J=1,N 
XBrtkl J)=T( J) 

150  T(J)=T(  JI/bK.K 
C 

C  CALCOLATb  SONS  Ob  Ck U SS- PkOOOC T S  Ob  DEVIATIONS 

C  bKlli-1  TEwRORAkY  NEAR'S  EUR  I'i  ( > BS bk V A T 1  OnS 

C 

L=U 

DO  lfio  1=1,  N. 

JK  =0 

mi  17u  J=l,k 
L  =  L  + 1 

17(1  I » (  J  )  =K  X  (  L  )  —  T  (  J  ) 

DO  180  J  =  1  , 

B  (.l)=o(  J)rD(  J) 

DO  1*0  K  =  1  , J 
J  K  =  J  K+  1 

lhO  k(  JK)=K(  JK)+D(  JX'D(K) 

C 

I  b  ( tv  —  K .<  )  203  ,  20  3,  18  3 
C 

C  kb  A  0  THb  kb  ST  Ob  OrtSbk  VA  I  Jlib  s  Oi\lf  AT  A  TINE,  SON 

C  THE  UnSEk  VAT  1 1  hv  ,  Aim-  CALCOLATb  SUNS  Ob  CkOSS- 

C  BkODUCTS  ob  DEVIATIONS  FkiJO  TbbPOkAkY  nEaOS 

C 

18  3  K  K  =  bi — K  K 

DO  200  1=1, KK 
JK  =0 

CAI.  L  DATA  (i",0) 

DO  190  J  =  1  ,  - 

Xb  AK  (  J  ) =X  h Ak ( J ) +0 ( J ) 

D( J)=D( J ) - T ( J) 

190, H( J)=b( J)+D( J) 

HO  200  J=1,D 
DO  200  K  =  1 ,  J 
JK  =  JK  +  1 

200  « ( JK )=k ( JK ) +D( J  )  *11  (  K ) 

C 

C  CALCULATE  n  E  A  n S 


COkkbl 1 2 
CORK  El 1 3 
C0kKF114 
C0RRF115 
C0RRE1 16 
CORK  FI  17 
CUkkbllB 
CORR  cl  19 
CO  Rk  b  1  20 
C0RRF121 
CO  k  k  E 1  2  2 
CORK  FI  23 
CO  kk  E 1 24 
CO  kk  E 1 2  5 
CUKkE l 26 
CO kk  F 127 
CO kk El  28 
CO  Kk  F 1 29 
CUkkbl 30 
CORK  F 1 31 
CORK  El  32 
CUkkbl 33 
CO  Kk  b 1 34 
CORK  b 1 3  5 
Ctl  k  K  b  1  36 
CORK  FI  3 7 
C0RKE13H 
CORK bl 39 
CO  K  K  b  1  40 
CORK b 141 
CORK  b 1 42 
CUkkbl 43 
CO kk  b 1 44 
CI)kkbl>+4 
C( 1  k  k  b  1 4  6 
CORK,  b  1 4 7 
Cork  b 1 48 
CUkkbl 49 
CUkkbl 30 
CUkkbl 31 
CUkkbl 32 
CO  KK  c 1 3  3 
(J  ■  kk  F  1  34 
CORK bl 33 
CUkkbl 36 
CUkkbl 37 
Cl  I  kk  b  1  38 
(  '  J 8 k  *  1 39 
COkkb 160 
CORK b 161 
CUKKE 162 

C  I  k  k  b  1  6  3 
Cllkkb  1  64 
COkkb  1  6  1 
( .  1 1 K  k  8  1  6  6 
Cl  IkK  b  1  6, 7 
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C 

201?  JK  =0 

no  210  j=i,p 

XBAR  (  J)=Xb*M  J)/Fi'J 
C 

ADJUST  SUMS  of-  CkUSS-PROi  JCTS  UF  DEVIATIONS 
FROM  TEMPORARY  MEANS 

DO  210  K=1,J 
JK  =  J  K+ 1 

2  10  K(  JK)=k(  .  ;<)-h(  J)*b(K)/FiM 

CALCULATE  COkk  F 1.  AT  1  ON  CoEFFl  S,I  FimTS 
JK  =0 

1)1)  22()  J  i,m 
JK= JK+ J 

220  STUio)=  SOKTI  ABS(K( JK) ) ) 

DO  230  J=1,m 
n r*  2  3o  k  =  j,i*i 
JK  =  .;+  i  K-K-K  )  /2 
[  =,.t*(  J-l  )+K 
KX(L)=R(JK) 

L.  ='•)*<  K-l»+J 
kX  (  L  )  :-K  (  JK  ) 

2^0  R  ( JK  )  =  «  i  JK  )  /  (  S TO (  ,. )  #S Tl) (  K  ) ) 

calculate  standard  deviations 

Fw  =  SDK  T  (  Fm-  1 . 0  ) 

D( i  2^0  ^  -  1  *  M 
240  s  T  0  (  J  >  =  STD  (  J  >  /Fix. 

C 

C  COPY  Tut  01  AOOXiAL  OF  THF  MATRIX  OF  SUMS  UF  CkUSS-PkOuUCTS  OF 

C  OF  V  i  A  T  J I II')  S  FROM  Mt«NS. 

c 

L  =  -m 

no  2  50  I  - 1 » t" 

L  =  L  +i-i  +  1 

250  B  ( I  )=«X (L ) 

KPT  UK  iv 
F  Ml  i 


CORK  FI  68 
CO RR FI  69 
CORK  FI  70 
CORK  FI  71 
CORK  F 172 
CORK  FI  73 
C0RRF17A 
CO RK  F 1  75 
CO  R  K  F  1  7  6 
CORK  FI  77 
CORK  FI  78 
CORK  FI  79 
CORK  F 180 
CORK  FIR  1 
CO  RR  F 182 
CORK F 183 
C0KKF184 
CORK  F  1  8 5 
CORK  F  1  HF> 
CORK F  187 
CORK  F 188 
Cl  IKK  F  1  8  9 
CflKKF  190 
CORK F  191 
CO  R  K  F  1 9  2 
CORK  FI 93 
CORK  FI  9  4 
CORK F) 95 
CORK)- 198 
CORK  F  197 
Cl  i  k  K  F  1 9  8 
C 1 1  k  K  F  1 9  9 
C0RKF200 
CORK  F?C  1 
Cl  I  PR  h-  ?()  2 
CURKF203 
CORK  F 20  A 
CO  RK  F 20  5 
CORKF70F 
Co  RK  F  20  7 
Cl  IF  KF  20  8 
CORK  F  20  9 
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SUBROUTINE  GP 


SUBROUTINE  GP  (X,  Y,  L,  S,  N,  I'l,  H,  LN,  A,  PLOT) 
CONTROL 


CaLL  UP  (X  ,  Y  ,  L  ,  S ,  R* »  W ,  L  N »  A  f  PLUT) 

lM  H  S  K  t- 

X  =  AkkAY  UP  I  NDEPENDFMT  VALUES,  U I D '«S  Iuupn  x(w). 

Y  =  AkkAY  UP  SFTS  UP  DEPENDENT  VALUPS ,  l>  ImEnS  Iuotl)  n(m,m). 
L  =  PIliHbEK  HP  LINES  ‘  i  6 1  SKiPPEI)  P  FU  kp  niSPLAY. 

S  =  iMlii-BPk  UP  SPACES  FROM  LET  SIU  i:P  PACE  TU 
BP  SKIPPPU  BEFORE  DISPL  Y. 

i-i  =  dunbek  points  iim  pach  set. 

i'l  =  NUMBER  UP  SPTS  OF  PU  IMS. 
w  =  WIDTH  UP  DISPLAY  IN  PKI1MT  Si-  ,ES. 

L  ■'A  =  L  EiMGTH  l)F  DISPLAY  IN  PRINT  L  .,vf  . 

A  =  AhKAY  UP  SINGLE  CHARACTERS,  DU'  N  S  I  ( ! P  U  A(n),TO 

k  EPKESEnT  Trip  TREND  FUR  EACH  St  (  F  X  .  —  UaTh  a/1hA, 
1HB,,.,ETC. ) 

PLUT  =  AkkAY  UP  SINGLE  CHARACTERS  GENERATED  BY  GP  TU 
DISPLAY  TRENDS,  DIMENSIONED  PLUT  iLN,W). 

Integer  s ,  w »  v-*  1 

I!  I  'JfcNS  I  UN  X  ( i'i )  ,  Y|H,N),  A  (  N  )  ,  PLUT(LN,W) 

DAT  A  PL  An  R  /  1 H  /,  pl'GP/lHf/ 

CHECK  MAX  I  WIDTH  Aim!)  LENGTH  kEDUFSTPU  and 
EXIT  IE  NI.lT  CURkbCT 


Ik,  p  LY  )i)  1 
DSPLYOO? 
i  iSp  L  you  i 
1)SPLV‘  )A 
1 1 S  P  L  ;  .  R 
I S  P  L  Y  0  0  f ' 
DSPLYOi)  7 

DSP  LYD1’1' 
I'SPlYI  id 
DSP  L  y()  10 
1 1 S  u L  V  11 
DSP  LYu  i  y 
DS--LYH  1 
DSPLYU  1  A 
DSP  L  YO  1  w 
I ‘S  P  L  Y;  Ik 
uSP  i_Yo  i 
DSPLYOI- 
1 1  S  P  L  Y  i  ’  1  L< 
I 'SR  L  V  t !  1 0 
USPLYu  L i 
D  i  p  L  Y 1 1  2  i 
US  p  LYi  1 
i.'S  ■  LYu  p 
US  LYU  <t  k 
l  'S  LYO/h 
i  I  S  p  i  v  u  ?  7 
I  s  i-  v !  v  p 

!'<(•  Yll,-'u 

' 1  s  P  L  ■  v ' 
'-  i  t  Vi'  -l 


IF  (S-iu  .gT.  131)  (it i  TU  9  00 
IF  ( L  +  L  n  ,GT.  S>k  )  GU  TU  HOD 

F  I  .\  1 1  m  1 N  1 1-1  U  Akil)  m«X  I  HU'''  UP  X  Ai'ki  Y 


I  '  S  ,  DIM 

u  S 1 1  L  v 1  1  y  * 

IMP  v,  -  i  ■ 


XNAX=X ( I ) 

X  H  1 1';  =  X  (  1  ) 

'I!  I  10  I=2,i- 

IP  (X  (  I  )  .GT.  X 1,1  A  X  )  *i.  AX-X  (I  > 

IP  (XU)  .LT.  XMM  X  i-1 1  l-'=  X  (  I  ) 

Y  maX =  Y  <  1,  1  ) 

Y  i-i  1 1-  =  Y  (  1 ,  1  ) 

DU  DU  l  =  l,<- 

DU  ?<)  >=1,N 

IP  (  Y  (  I  ,  J  )  .  G  T  .  Y-AX)  Y  k-AX  =  Y  !  1  ,  J  ) 
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APPENDIX  IV 
STATISTICAL  FORMULAS 


Mean  (Arithmetic  Average)  -  M 

"i 

Mi  3  i 

i  :| 

where 


=  each  individual  value  (SN)  in  a  column  "j" 
=  total  number  of  X„  in  a  given  column  "j" 


Standard  Deviation  -  SD 


SD  is  a  measure  of  dispersion  ("scatter")  of  a  given  column  of  data. 


Coefficient  of  Variation  -  CV 

so 

cvi  =100  x  — — 

1  u 

"i 

Since  CV  is  a  calculation  of  dispersion  per  unit  mean,  it  is  an  excellent 
indicator  of  precision.  Minimization  of  CV  is  the  goal. 

Pri  duct -Moment  Correlation  Coefficient  -  r 

<7^  1  t  ud-“i  >  > 

i  •  I 

V  ( SO  j  )  (SOt  ) 

r  between  any  two  columns  of  data  "j"  and  "k,* 
n  *  n^  *  n^  =  to  al  number  of  data  points  per  column. 
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This  "r"  is  often  called  the  sample  correlation  coefficient.  This  can  be 
related  to  "  P, "  the  population  correlation  coefficient,  as  a  function  of  sample 
size  n  for  given  confidence  limits,  by  using  standard  graphs  (Reference  6). 

The  theory  and  derivation  of  these  quantities  can  be  obtained  from  most 
textbooks  on  engineering  statistics,  including  Reference  9. 
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APPENDIX  V 

IMPROVING  SMOKE  MEASUREMENT  RESULTS  TAKEN 
WITH  MILLIPO1  MEDIUM 

As  explained  in  Appendix  II,  the  smoke  measurement  (SN)  is  a  function  of 
Q,  the  sample  size  per  unit  filter  media  area  (scf/sq  in). 

All  other  factors  being  constant,  Q  is  proportional  to  the  amount  of 
particulate  matter  per  unit  filter  media  area,  "W"  (micrograms/sq  in).  Also, 
by  definition,  SN  is  a  function  of  the  absolute  reflectance  of  the  smoke  spot. 

As  also  explained  in  Appendix  II,  loading  curves  for  given  filter  medium 
are  graphical  relationships  of  spot  reflectance  (or  optical  density)  versus  W. 
Since  these  are  exactly  analogous  to  SN  vs  Q  functions,  this  author  spent 
several  hours  working  with  Millipore  and  Whatman  media  loading  curves  pro¬ 
duced  by  GE  in  an  attempt  to  find  a  way  of  improving  the  results  that  would  be 
obtained  with  Millipore  medio,  "Improving"  in  this  case  means  increasing  the 
average  SN  of  Millipore  to  the  same  level  produced  by  Whatman  medium  (the 
analysis  conclusively  demonstrated  that  Millipore  produced  results  consistently 
lower  than  those  of4  ained  with  Whatman  medium). 

The  SN  from  Whatman  and  Millipore  media  are  reported  at  Q  values  of 
0.  300  and  0.0565  scf/sq  in,  respectively.  Since  "N  varies  withQ,  the  problem 
was  to  find  that  vaiue  of  Q  for  Millipore  medium  that  produced  the  same  SN 
as  the  accepted  value  of  Q  for  Whatman  medium  (0.  300  scf/sq  in).  To  do  this, 
the  author  worked  with  average  differences  in  SN  taken  from  Table  I  and  the  GE 
loading  curves. 

Such  an  approach  can  yield  an  approximate  answer,  at  best,  because  the 
difference  between  the  two  loading  curves  is  not  constant  (one  curve  is  not 
merely  the  translation  of  the  other).  The  curves  diverge  increasingly  with 
increasing  W.  Figure  9  demonstrates  this. 

Only  the  lower  portion  of  the  SN  range  (about  10  to  50)  was  used  to  mini¬ 
mize  the  effect  of  this  divergence.  This  is  reasonable  since  SN  of  about  10  to 
50  are  of  greatest  concern  (see  Section  1112). 
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The  final  answer  was  that  Millipore  with  a  Q  of  0.  0660  scf/sq  in  would  give 
SN  of  about  the  same  as  Whatman  wit  =  0.  300  scf/sq  in.  It  cannot  be  over¬ 
emphasized  that  thi3  is  an  approximate  answer.  Additional  experimentation  is 
necessary  to  corroborate  this  value. 

Even  if  Q  =  0,  0650  scf/sq  in  proves  to  be  the  "proper”  value  for  use  with 
Millipore  medium,  it  will  still  be  accurate  only  for  part  of  the  range,  although 
fortunately  the  most  important  part  (SN  of  10  *o  50).  It  shouH  also  be  noted  that 
this  adjustment  of  Q  affects  the  magnitude  of  resultant  SN  only;  the  greater 
dispersion  of  Millipore  results  will  probably  still  be  present. 
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